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A surface ship subjected to an underwater explosion is exposed to shock waves 
over a short period of time which can vary in magnitude based on charge type, size, and 
location. The energy of those waves impinging upon the hull is transmitted throughout 
the ship’s structure and vital equipment. The dynamics of the shock waves also influence 
the fluid surrounding the outer hull of the ship, creating an area of cavitating fluid. The 
combination of the shock waves, bubble pulsations, and cavitating fluid induce shipwide 
vibrations on hull supports and mission essential equipment which may become 
inoperative. In view of congressional requirements for new ship designs and systems to 
be shock tested, this thesis investigates the modeling of a preliminary design (Flight I) of 
the Arleigh Burke Destroyer (DDG 51) exposed to an underwater explosion. The effects 
of cavitation on one and two dimensional models is explored to determine if cavitation 
effects are substantially important to a three dimensional ship model. Validation of 
modeling underwater explosion effects upon a ship model can provide potential insight 
and savings in cost for future live fire testing and evaluation of the Flight IIA (DDG 79) 
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1 . 


INTRODUCTION 


A. BACKGROUND 

Today’s U.S Navy consisting of technologically advanced ships such as the 
Ticonderoga class cruiser and the Arleigh Burke clciss destroyer are being driven to 
operate in the littoral, close to shore, environment as dilineated in the Department of the 
Navy’s strategic planning document [Ref. 1]. These platforms equipped with the Aegis 
weapon system were designed for open ocean, “blue water” tactical scenarios but due to 
the end of the cold war, a shift in emphasis is placed on low intensity conflicts which 
occur near the sea-land interface. This has increased the potential for “cheap kill” 
opportrmities against U.S. weirships such as an underwater mine as ships close the shore 
and the availability of such weapons is high and the cost low for third world country to 
purchase. 

As recently as 1991 during Operation Desert Storm, the USS Princeton (CG-59), a 
Ticonderoga class guided missile cruiser, struck a floating mine near the bow. Although 
the rupture to the hull was localized, the subsequent whipping of the hull caused severe 
hull girder damage near the stem of the ship and weapon systems to go off line. This 
type of tactical environment now places a greater responsibility on ship designers and 
engineers to meet requirements as specified in MIL-S-901D [Ref. 2]. 

A number of studies have been conducted on submerged cylinder type models 
exposed to an underwater explosion in order to predict the hull response of stmctures 


such as a submerged submarine. These studies did not include the effects of cavitation on 
the wet surface of the model or of the cavitation which occurs at the free surface. For 
surface ship type studies, simplified beam models such as that conducted by Hicks [Ref. 
3] have concentrated only on the effects of whipping of surface vessels due to gas 
bubble dynamics. In this age of improving computer hardware, larger memory and 
software programs such as the USA (Underwater Shock Analysis) code with the CFA 
(Cavitation Fluid Analyzer) [Ref. 4] to model underwater shock effects and improving 
CAD programs with built in FEM (Finite Element Model) generators such as IDEAS 
[Ref. 5] and MSC/NASTRAN [Ref. 6], the engineer has a greater ability to improve ship 
modeling and predict the vibrational response of ship hulls and internal equipment. 

B. SCOPE OF RESEARCH 

This thesis studies the methods of modeling one and two dimensional finite element 
models with and without an external fluid. Two different computational methods are 
selected in the USA code for each model arrangement in addition to incorporating the 
CFA module in the solution process. The goal is to determine the cavitation effects on a 
surface vessel and develop the most simplest and conservative method in modeling the 
underwater explosion. A one dimensional model is created from Bleich’s paper [Ref. 7] 
which defines the conditions of a square plate sitting on a fluid surface and a fluid 
column. In additon, an exponetial decaying plane wave models the underwater explosion. 
A second finite element model is created to represent the cross-sectional area of a simple 
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ship hull. This two dimensional model includes the surrounding fluid to the waterline of 
the hull and introduces the effects of a free surface due to an actual charge placed at a 
distance from the hull. The results from these first two models provides some 
simplification for the three dimensional notional hull structure of the Flight I of the 
Arleigh Burke. Although live fire testing and evaluation has already been conducted on a 
Flight 1 design, a Flight IIA design starting with DDG-79 possesses hull differences 
shown in Fig. 1 which may require shock testing. Improvements in the area of modeling a 
surface ship such as the Flight IIA design from existing CAD programs and fluid 
behavior during an underwater explosion through the USA code, can create potential 
savings in future Flight IIA shock trial tests. 



Figure 1. Major Changes (Flight II to Flight IIA) [Ref. 8] 
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11. UNDERWATER EXPLOSION 


A. FLUID BEHAVIOR 

When an underwater explosion occurs, the solid explosive material (such as 
HBX-1, TNT, etc.) suddenly reacts, leaving behind high temperature and pressure of 
gaseous products. An initial wave of compression is produced due to a sudden 
discontinuity of pressure which can reach a peak pressure of 2 x 10^ psi. This steep 
fronted wave known as a “shock wave” propagates radially outward at a speed several 
times the speed of sound in water and roughly decays exponentially. The velocity of the 
shock wave decreases such that at approximately 10 times the explosive charge radius, 
the disturbance is essentially traveling at the acoustic velocity in water, 5000 ft/sec. [Ref 
9:pp 3-7 ] 

The pressure sensed at a specific point in the fluid or at a structure is essentially a 
step increase to the peak pressure followed by an exponential decay with time until the 
pressure has dropped to approximately one third of the peak pressure. After falling to 
approximately one-third of the peak pressure, the pressure level falls off inversely with 
distance [Ref 9:pp3-7]. 

The cavity of gaseous products left behind at high pressure forms a bubble which 
subsequently expands to relieve the difference in pressure, accelerating the surrounding 
fluid particles. The bubble continues to expand beyond the point of hydrostatic 
equilibrium (due to the inertia of the surrounding fluid) until a point of dynamic 
equilibrium is reached. The bubble then reverses its motion, continuing to contract until 
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dynamic equilibrium is again reached, where it quickly rebounds and again begins to 
expand as illustrated in Fig. 2. This oscillating bubble expansion and contraction 
continues until the energy of the reaction is fully dissipated or the bubble finally reaches 
the surface, venting the gaseos by products of the explosion. As the bubble rebounds, it 
greatly accelerates the surrounding water, generating a substantial pressure pulse (known 
as the bubble pulse). This bubble pulse can impart significant loads on structures in the 
vicinity. 

In addition to the initial shock wave and bubble pulse, based on the location of the 
expolsive charge with respect to the surface and sea floor, a surface vessel can also 
experience a combination of the following pressure waves as illustrated in Fig. 3: 

1. Free-surface reflection 

2. Bottom reflection 

3. Bottom refraction (not shown) 

The directions of free-surface and bottom reflection waves can be determined through an 
application of Snell’s law at the surface and sea bottom, respectively. At the surface of a 
fluid, a reflected wave of negative pressure known as the “rareffaction wave” is formed 
with a value such that the sum of the direct and reflected pressures is zero along the 
boundary condition between air and water defined by the surface. This rareffaction wave 
travels through the fluid region shortly after the incident shock wave at any point which 
causes the incident shock wave pressure profile to be truncated at a point in time called 
“surface cut-off.” Bottom reflected and refracted waves are dependent on the 
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characteristics of the sea bottom and can contribute significant pressure waves in shallow 
water environments but are not considered for this study. 



Figure 2. Gas Bubble behavior 



Figure 3. Pressure Wave Profiles 
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Figure 4. Combined Pressure Profile of an Underwater Explosion 

The oscillating bubble pulse can also introduce a low-frequency forcing function that can 
resonate a ship girder frequency especially when the period of the bubble pulse is in the 
vicinity of the natural frequencies of the ship girder. Fig. 4 illustrates a combination of 
the pressure waves. 

A number of empirical relationships have been derived through numerous 
underwater explosion tests since World War II which provide a means to calculate 
explosion characteristics such as incident wave pressure (Pmax) which impinges upon a 
target, shock wave decay constant (0), pressure as a function of time (P(t)), and bubble 
period (T) with derived constants for each respective explosive material provided: 
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W = Charge weight (lb) 

R = Standoff distance (ft) 

D = Charge Depth (ft) 
ti = Arrival time of shock wave (sec) 
t = Time of interest (sec) 
K|,K 2 ,K 5 ,Ai,A 2 = Shock wave parameters 


B. STRUCTURAL BEHAVIOR 

The dynamic response of a linear-elastic structure can be defined, based on the 
number of elements and nodes that define the structure, as 

[M]{x(t)}+[C]{x(t)}+[K]{x(t)}= {F(t)) (2.1) 
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where [M] = structural mass matrix, [C]=damping matrix, [K]=stiffiiess matrix, 
{F(t)}=applied external force, {x(t)}= general displacement vector and derivatives with 
respect to time. For a structure in a fluid medium, the applied external force created by an 
underwater explosion, is defined as 

-[G ] ({Pi(t)} + {Ps(t)})= {F(t)} (2.2) 

where [G] = the transformation matrix between structural and fluid nodes, [Af]ui(i]= 
diagonal area matrix for the elements in the fluid, {Pi(t)}= incident wave nodal pressure 
vector and {Ps(t)} = scattered wave nodal pressure vector. 

The USA code provides the means to model the pressures involved in an 
underwater shock as a “forcing function” on a fully or partially submerged model. The 
code also utilizes variations of the DAA (Doubly Asymptotic Approximation) method 
[Ref. 10] that models the surrounding acoustic medium as a membrane on the surface of 
the structure actually in contact with the homogeneous fluid. The effects of cavitation on 
a structure modeled with a surrovuiding fluid can be also included through the CFA. The 
fluid motion and structural response is linked together through a series of compatibility 
relations that are matrix differential equations that relate the symmetric fluid mass matrix 
of the model’s wet siuface, fluid particle velocities from the scattered wave normal to the 
structures wet surface, and rates of changing scattered fluid pressures. The USA code 
creates this matrix by utilizing an infinite, inviscid, incompressible fluid in an irrotational 
flow created by the structure’s wet surface and then creates a virtual mass matrix which is 
comprised of the added mass of the fluid on the elements of the structures mass matrix 
which are wet. Kinematic compatibility is obtained between the fluid motion and 
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structural motion by an invariance of virtual work and transformation matrices that 
create the final interaction equations for fluid and structure motion [Ref 10]. 

Hence, the final equations of structure and changing fluid pressures create a set of 
“ augmented interaction equations” which the USA code solves, given a specified time 
increment, in an unconditionally stable manner. 


C. CAVITATION 


Two types of cavitation which can occur during an underwater explosion is that of 
“local cavitation” which occurs at the fluid-surface interface and “bulk cavitation” which 


occurs near the free surface that can extend beyond the target of interest by hundreds of 


feet. Utilizing Taylor flat plate theory, the particular models in this thesis will utilize only 


air backed plate conditions as illustrated in Fig. 5. 


A i r 
Backed 




Pi(t) 


Vi(t) 


Vp(t) 


Air 

Backed 



Pr(t) 


Vr(t) 



Shock Wave Approaching Shock Wave Reflecting 


Figure 5. Taylor Plate Subjected to Plane Wave. 
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1 . 


Local Cavitation 


When the pressure pulse from an underwater explosion impinges upon a flexible 
surface such as the hull of surface ship, a fluid-structure interaction occurs which 
dynamically excites the ship structure. Looking at an infinite, air backed plate of mass, 
per unit area (m) subjected to an incident plane shock wave Pi(t), a reflection wave of 
pressure P/t) will come off the plate. Letting Vp(t) be the velocity of the plate and 
applying Newton’s second law of motion: 

mx-^ = P, + P, (2.3) 

Letting the fluid particle velocities behind the incident and reflected shock waves be Vi(t) 
and Vr(t), respectively, the velocity of the plate becomes 

Vp(t) = Vi(t) - v,(t) (2.4) 

Incident and reflected shock wave pressures are defined, respectively : 

PrpCvi (2.5) 

Pr= pCVr (2.6) 

where p=fluid density and C=acoustic velocity. Substituting the pressures (2.5) and(2.6) 
into (2.4) and utilizing for a high explosive shock wave 

Pi(t) = P™.xe® (2.7) 
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where t = time starting after the shock wave arrives at the target, one can solve for P^: 


P, = Pi -pCVp = ® - pCVp 

Returning to the basic equation of motion (eq. 2.1) and substituting for Pj gives the 
following first order linear differential equation: 


dv 

^ -1" pCv„ = 

dt 


(2.9) 


The above differential equation is solved for Vp, 




2P.axe 


m(l-p) 


e ^ —e ^ 


( 2 . 10 ) 


where P=pC0/m and t>0. Total pressure on the plate is also found to be 


P. + P = P 

I r ^ max 


2 A 2P At 


•e ® - 


1-p 1-P 


( 2 . 11 ) 


Hence, as (3 becomes larger, representing a light weight plate becoming more 
flexible, the total pressure will become negative at a very early time and as the total 
pressure approaches vapor pressure, local cavitation occurs in front of the plate. At this 
point, the plate is separating from the fluid and the maximum velocity for the plate is 
attained [Ref. 11] 
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2. 


Bulk Cavitation 


The combination of the rarefraction wave and incident shock wave near the free 
surface creates a region of decreasing fluid pressure that approaches the vapor pressure of 
water. The rarefraction wave direction is determined by an image charge located an 
equivalent distance from the surface boundary in the air as the actual charge in the fluid. 
An upper and lower boundary begins to form which defines the bulk cavitation zone (Fig. 
6). The velocity of the water particle near the free surface will have primarily a vertical 
velocity at cavitation. As the reflected shock wave passes, the fluid is then acted upon by 
gravity and atmospheric pressure. 



Figure 6. Bulk Cavitation Zone [Ref. 11) 
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The upper cavitation boundary of Fig. 6 is the set of location points where the rarejfraction 


wave arrives and the absolute pressure drops to the cavitation pressure. The total pressure 


which determines the upper boundary is the summation atmospheric, hydrostatic, incident 


shock, and rareffaction which can be defined as 


F(x,y) = K,(-)^‘e C6 +p^+yy_K,( 

r, 





r^=^{D-yf +x^ 




where 


x,y = horizontal range and vertical depth of point 


ri = standoff distance from charge to point 
Xj = standoff distance fi'om image charge to point 
C = speed of sound in water 
D = charge depth 
0 = decay constant 
Pa = atmospheric pressure 
y = weight density of water 
W = charge weight 


Ki,A, = shock wave parameters 


The lower cavitation boundary is determined by equating the decay rate of the 


breaking pressure which is the rarefaction wave pressure that reduces a particular location 


of fluid to the point of cavitation pressure, approximately zero, to the decay rate of the 
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total absolute pressure of the fluid as illustrated in Fig. 6. This boundary can be defined 
as [Ref 12] 


G(x,y) = -^ 
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R J 
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Ai 


Hence, any point that satisfies F(x,y)=0 and G(x,y)=0 determines the upper and 
lower boundary points. Fig.7 and Appendix A provides a MATLAB code [Ref. 13] for 
bulk cavitation zones for various depths and charge weights for TNT explosive material. 
This particular tool was useful in modeling the 2D cross sectional hull and placing the 
charge in location that would minimize the effects of bulk cavitation and aid in 
concentrating effects of local cavitation. Fig. 7 illustrates the region of the bulk 
cavitation zone above the charge does not vary much with changing depth. However, the 
bulk cavitation zone increases in depth with appreciable increases in charge size. 30 feet 
below the keel (total depth of 50 feet) was used for the two dimensional study. 
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Depth(FT) Depth(FT) Depth(FT) 


100LBS TNT 





Figure 7. Bulk Cavitation Zones for a TNT charge located at the following depths 

—50ft,-100ft, -.150ft 
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III. COMPUTATIONAL FLOW 


Model preparation, translation, calculation, and post-processing involved a 
number of programs which is depicted in Fig. 8. Each process involves a number of 
output files which are input to the next step in the process or provides the user a means to 
verify results from each step before proceeding. 


MODELING 


TRANSLATION 


USA SOLUTION 


TRANSLATION 


POST-PROCESS 



Figure 8. Computer Program Flow 
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A. 


MODELING 


1. ASSET 

ASSET/MONOSC (Advanced Surface Ship Evaluation Tool/Monohull Surface 
Combatant) [Ref. 14] program is an interactive program which created the original ship 
dimensions for the three dimensional study. It provides a naval architectural tool to 
address ship design, which includes geometric definition of the hull and superstructure, 
hull subdivision, hull structure, resistance, appendages, equipment, and costs [Ref. 14]. 
An early design of the DDG 51 class destroyer with general weights of equipment 
groups and hull structure requirements based on inputs from designer imposed needs or 
mission requirements was exported out of ASSET in an IGES (Initial Graphics Exchange 
System) file format which included general ship offsets which can be imported into 
IDEAS through the IGES translator. Although general bulkhead, deck, superstructure 
design was provided, no further details to subcompartmentized spaces, specific weights 
and locations of internal equipment, sonar dome offsets, and any other explicit details of 
the notional ship design was available. Appendix B provides a summary of the DDG 51 
notional model. 

2. IDEAS 

IDEAS, Master Series 2.0, [Ref. 5] was utilized to model all three cases of study. 
The IDEAS program is comprised of a large number of applications such as design. 


20 


modeling, testing, manufacturing, and simulation, each with a large group of subtasks for 
the user to utilize in icon panel form or menu driven selection. The program is also a 
“team concept” software which allows a number of users in a design or analysis team to 
take individual parts of the whole model to be designed, modified, or analyzed. For each 
model, the Simulation Application and all its respective subtasks were utilized to create 
(e.g. surfacing) and mesh the finite element model for export. 

a. One and Two Dimensional Models 

For the one dimensional and two dimensional cavitation study, one model 
was required to be created which included the structure and surrounding fluid. Fluid mesh 
thickness could not exceed the following requirement due to numerical stability of the 
CFA code; 

iP^<5 

P.'. 

where =density of water, pj =density of steel, tj =thickness of steel plate, and D= width 
of fluid elements. For example, utilizing p/p^ « 8, if the thickness of the plate or hull % 
=0.5 in, D < 10 in. After creating all structure and fluid model regions, the following 
order for meshing was conducted due to the order that the CFATRAN and CFAGAL 
program organizes particular fluid nodes and elements: 

1. Fluid volume elements (brick elements) 

2. DAA boundary (thin shell elements) 

3. Hydrostatic boundary (thin shell elements) 
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4. Specified Pressure regions (not utilized for this study, thin shell elements) 

5. Wet structural surfaces (thin shell elements) 

6. Remaining structural surfaces 

The thin shell elements for the DAA boundary, hydrostatic boundary, as well as for the 
specified pressure regions if used, provide a means to select the nodes in these particular 
regions that the CFA requires for numerical solution. The thin shell elements do not 
become a part of the soluiton process but act as “markers” in the CFATRAN program for 
the boundary element nodes which are specifically associated by the MAT field in the 
NASTRAN ASCII input file. Because of the numerical processing in the USA code, 
node numbering needs to be described in a counter-clockwise manner looking from the 
fluid towards the structure. Color coding particular surfaces and their respective elements 
or by creating specific groups of surfaces or elements in the meshing subtask-grouping 
icon greatly aids in changing particular characteristics such as the node numbering. 

A one psi pressure vector outward from the structure into the fluid is 
applied to all structural wet elements to act as a marker for the USA code fluid matrix 
development and restraint set conditions set for symmetric boundaries are created in the 
BOUNDARY CONDITION subtask. 

After meshing, two copies of the global mesh are made. One FEM model 
will be structiue only with the one psi pressure, deleting fluid elements and nodes and 
one FEM model will be the fluid volume elements with the DAA, hydrostatic, and wet 
structure thin shell elements, deleting the other dry structural nodes and elements. These 
two FEM models are then exported out of IDEAS and import to MSC/NASTRAN. 
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b. 


Three Dimensional Ship Model 


Prior to the IGES translator task, modifications needed to be made to the 
actual IGES ASCII file for the correct unit system, ensuring IDEAS was preset to the 
same units (feet, inches, meters, etc.). 

Offset points and lines .were utilized as “wire guides” to recreate the model 
in IDEAS for creating all surfaces through lofting, mesh of curves, or planar creation 
shown in Fig. 9. The following sequence proved to be the best option in creating a 
complete three dimensional ship model in IDEAS, Simulation, Master surfacing task: 

1. HULL 

2. BULKHEADS 

3. DECKS 

4. SUPERSTRUCTURE: 

After all the surfaces have been created, all the edges must be connected or 
STITCHED in order for meshing to be successful and a check for free edges is performed 
in preparation for meshing. Any new surfaces that need to be created, such as internal 
bulkheads or connecting the superstructure to the hull, can be done through a surface 
bound by a combination of edges and three dimensional lines connecting those opposing 
edges. 

The required thin shell element properties are applied and a one psi 
pressure vector must be applied to all the wet surfaces through the BOUNDARY 
CONDITION subtask as well as any symmetry condition through the RESTRAINT SET 
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icon. IDEAS provides a full series of element checks ('warping, distortion, interior angles, 
and normal directions) which must be done for proper solution processing. Only one 
FEM model was required for export since the CFA program was not applied to this 
particular model. 
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Figure 9. Imported ASSET DDG 51 Model 
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B. TRANSLATION 


Each of the following steps allows for output files to be viewed, verifying that 
information (number of nodes, eigenvalues, etc.) is being properly computed and created 
in the unformatted fortran files for the next step. 

1. MSC/NASTRAN 67.5 

The MSC/NASTRAN structure database is the required format for this particular 
version of the USA code during the initial start and restart. The USA code itself is 
conducting the solution process which would normally be a SOL 109 transient analysis 
process in MSC/NASTRAN. Appendix B provides the DMAP alters for IDEAS post 
processing which is attached to the beginning of the exported file out of IDEAS and for 
the MSC/NASTRAN restart, deleting unnecessary and unneeded lines up the bulk data 
section. This process generates the geometry information and the structure’s M, C, and K 
matrices 


2. CFATRAN/CFAGAL 

CFATRAN is a program modified from [Ref. 15] that takes the 
MSC/NASTRAN file of the fluid mesh model and groups the elements, nodes, and nodal 
connectivity in the required format for the CFAGAL program. CFAGAL will read the 
CFATRAN output and create a global access library of the fluid model and its 
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boundaries, just as the program FORGAL does for the structure model. This library will 
then be accessed during the FLUMAS process. 

3. NASSKY/FORGAL 

Nassky prepares the NAS TRAN UTl output file which contains the geometry and 
matrix data of the structure into “skyline” matrix format, a generic Fortran unformatted 
file. FORGAL then converts this generic file data into a GAL file (Global Access 
Library) database which the USA code will access for structure related data. 

C. UNDERWATER SHOCK APPLICATION AND SOLUTION 

Appendix B provides sample inputs for the following programs which were 
utilized for the two-dimensional study. 

1. FLUMAS 

This is the first step in the USA code which creates the fluid mass mastrix of a 
structure wet surface in an infinite, inviscid, incompressible fluid by the boundary 
element technique. The key inputs here are the following [Ref. 16,17]: 

1. Mesh Geometry: fluid wet-surface and structure wet-surface mesh 

2. Element definitions: general curved surface or surface of revolution 

3. Material Property: mass density, speed of sound, DAA2 parameter 

4. Constraint options: location of free surface, half model, quarter model, long 
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cylinder, node reassignment in fluid-structure transformation if needed. 

For cavitation studies, the CFA is activated and the fluid volume database created 
by CFAGAL is incorporated with the FLUMAS output. Care should be exercised here 
during CFA solutions. The NCAV field refers to the number of elements on the wet 
structure plus the DAA boundary. During non-CFA runs, the NCAV field refers only to 
the wet elements on the structure since this defines the DAA boimdary. 

2. AUGMAT 

This process in the USA code utilizes the equations specified earlier as the 
“augmented equations” by utilizing the data in the GAL file and FLUMAS output to 
construct the specific matrices. The key inputs are the following [Ref. 16,18]. 

1. Fluid and Structure mass matrices 

2. Structural DOF correspondance table 

3. Fluid mesh geometry 

4. Fluid/Structure DOF transformation coefficients 

5. Type of solution to be conducted; plane wave, DAA2 options 

3. TIMINT 

This step constitutes the most time consuming step depending upon on the 
number of nodes, elements, and requested length of time to be studied. The 
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unconditionally stable staggered soluiton technique is utilized which basically involves 
the following [Ref. 16,19]: 

1. Estimating the unknown structural restoring force [K^] {x(t)} at t + A t from the 
extrapolation of current and passed values. 

2. Solve the fluid equation and obtain a preliminary estimate of the total pressure 
vector at t + A t. 

3. Solve the structural equation for the displacement and velocity vectors at t + A t. 
Other steps are involved between these particular steps for matrix transformations and 
improved accuracy. General inputs from the user involve the following: 

1. Incident wave characteristics: location and type of charge source, standoff, 
specified pressure profile, exponentially decaying wave. 

2. Time step information: start and finish times, increments. 

3. Display data of particular nodes for checking output results in the solution 
file:desired node number and freedoms for displacement, velocities, or 
pressures on the DAA boundary. 

Although unconditionally stable, the user should utilize small enough time steps to 
capture accurate model behavior. Solution time history responses of all structural nodes 
and fluid nodes for all runs are stored in the file designated POSNAM in the TIMINT 
input. 


29 


D. POSTPROCESSING 


Two methods of post processing was done due to the growing amounts of 
required memory for larger models such as the three dimensional ship hull. 

1. XUDVT 

This second interface program converts the USA solution history in the 
corresponding MSC/NASTRAN data block, UDVT, a file format readable by 
MSC/NASTRAN. The initial MSC/NASTRAN job is restarted at the point just before the 
NAS SKY translation, allowing MSC/NASTRAN to continue to a normal conclusion. The 
USA solution replaces that which would normally be generated by MSC/NASTRAN in a 
transient analysis. 

During CFA studies, the number of wet points required to be entered are those 
wet points totals generated by NASSKY which can be found in the nassky.out file. 
However, during non-CFA studies, the wet points refer to the total number of wet 
element faces. If the incorrect number of wet points or wet elements are entered, 
translation will fail or data history for some nodes will not be translated. XUDVT also 
allows displacement, velocity, acceleration, and pressure histories to be extracted in 
ASCII table format and placed in the PLT file for X-Y plotting by other programs (e.g. 
MATLAB). 
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2 . 


POSTPR 


This is a USA postprocessing program which allows extraction of particular node 
displacement, velocity, or pressure histories. Although not needed for non-CFA studies, 
this is the only current method which one can extract pressure history profile of fluid 
nodes during a CFA study. The data is placed in ASCII column tables which can be 
plotted by other programs such as Matlab. Due to the large solution files generated by the 
three dimensional ship model and the growing size of files during the XUDVT- 
MSC/NASTRAN-IDEAS translation and final post-processing, POSTPR proved to be a 
memory saving tool. The user must ensme the original UTl file name is either deleted or 
changed in order to prevent a XUDVT failure since this process produces a UTl file for 
the MSC/NASTRAN restart. 

3. MSC/NASTRAN (restart) 

Appendix B illustrates the required fields needed to restart MSC/NASTRAN with 
a new UTl file created by the XUDVT process. This process utlizes the original 
model.MASTER file originally created by the initial MSC/NASTRAN run and creates a 
solution file for IDEAS to import (model.OP2) which contains all or selected time 
histories of displacement, velocity,acceleration, and stress. The user should utilize a 
different .OP2 file name from that of the .OP2 file created in the first MSC/NASTRAN 
run. Files usually triple in size compared to the POSNAM file output from TIMINT. 
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4. 


IDEAS 


The final .0P2 file can be imported into IDEAS through the FILE IMPORT- 
MSC/NASTRAN process. IDEAS will create a universal file (.unv) prior to actually 
creating a model file database. This unversal file is usually three times the memory file 
size of the .OP2 and model.POS file due to being in ASCII format. Model file database 
size (file.sfl and file.sfZ) become approximately the same size as the universal file. 
Hence, a user should ensure plenty of hard drive memory is available during this 
translation phase, especially for large models (> 2000 elements). 


32 


IV. ONE-DIMENSIONAL ANALYSIS 


A. PLATE MODEL 


A one dimensional study was conducted utilizing the case defined by [Ref. 7] 
which involved modeling a floating plate on the surface of a bilinear fluid. The goal of 
this model was to accomplish the same results shown in [Ref. 7] by creating a finite 
element model and solving with the USA code. Follow-on study was performed by 
Newton [Ref 20] in area of cavitation in support of developing a “shock shield,” a gas- 
filled cushion that was to surround the exterior of a submarine hull and reduce the 
magnitude of imderwater shock loads transmitted through the hull. 

Fig. 10 and Table 1 describe the particulars of the Bleich-Sandler case. 


1 thin shell element, 4 nodes, 1 ” thick 
Symmetry on all four edges 
Dimensions: 1.5”x 1.5” 

E=30.0 X 10^ psi 

Poisson’s ratio=0.3 

Mass density=5.32986 x 10"^ Ib/in^ 

Table 1. Plate Characteristics 


Fluid Depth: 150.0 in 

100 volume elements, 8 nodes, 1.5”X1.5” 1.5” 
404 total fluid nodes 
Mass density= 9.3455 x 10’^ lb 
Speed of sound= 5.712 x lO'* in/sec 

Table 2. Fluid Characteristics 
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Figure 10. IDEAS Bleich Plate and Fluid 
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Peak pressure =103 psi 
Decay time= 9.958 E-4 sec 

Charge depth = 1.0 E7 in, centered on the center of the plate 

Table 3. Incident Pressure Characteristics 


B. RESULTS 

Utilizing a time step size of 1.313 x 10'^ sec, 1200 time steps were collected. 
With the CFA interface turned on with the USA, program Fig. 11 depicts the computed 
velocity and pressure at the plate. Fig. 12 depicts pressures at specific locations along 
the column of fluid. Because of the numerical computation scheme within the USA code, 
time zero noted on each pressure figure refers to the point in time that the incident 
pressure wave arrives at that particular fluid node and not the global time of plate motion. 
Total pressures at time zero for each fluid node is the sum of the atmospheric, 
hydrostatic, and incident pressure wave at that point of arrival from the initial charge 
location and applied exponential decay. Fig. 11 illustrates that the cavitation region does 
not touch the interface boundary between the fluid and the plate. Node 401 in Fig 12 
confirms the cavitation region is captured by the fluid FEM and does not extend beyond 
the DAA boundary. Fig. 13 depicts the plate velocity and pressure with CFA OFF. 

Calculated time periods of zero pressure or cavitation occurrence grow to a period 
of .008 sec half-way through the column and stay approximately the same to the end of 
the fluid column. In addition, the CFA captured the cascading effects of cavitation 
closure which occured at similar times following incident pressure wave arrivals at each 
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fluid node and closure pressures became larger at deeper depths than at the surfaces. This 
closure phenomenon begins at the point t=.0108 sec at which the plate experiences an 
immediate velocity cutoff of the downward direction back to zero by the arrival of this 
“secondary shock” at the plate. Pressure at this particular time jumps to approximately 
14% of the peak pressure and creates an upward velocity of approximately 16% which 
brings the plate velocity back to zero. Multiplying time by 1004.22 to get it in terms of 
decay time units of the incident wave and normalizing velocity to 1/1000 of the speed of 
sound with a sign change due to the z-axis for verical displacement being defined 
downward into the fluid, calculated results match those (Fig. 14) with Bleich and Sandler 
as shown in Fig. 15 and 16. 

Hence, the finite element modeling approach with the USA/CFA code produced 
identical results as those shown by Bleich and Sandler through the method of 
characteristics [Ref 7]. In addition to the pressure history profile of Fig. 15, USA 
postprocessing was able to display the cavitation closure anywhere in the fluid column. 
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Figure 11. Bleich Plate Velocity and Pressure (CFA ON) 
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Figure 12. Bleich Fluid Pressure at Nodes 9, 200, and 401. (CFA ON) 
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Figure 13. Bleich Plate Velocity and Pressure (CFA OFF) 
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Normalized Vertical Velocity 



Figure 14. Normalized Vertical Velocity (—CFA ON, CFA OFF) 
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t, non-dimensional time 
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Figure 15. Time history of cavitated region (Ref. 7] 




Figure 16. Non-dimensional upward velocity ol plate [Ref. 7] 


41 













42 







V. 


TWO DIMENSIONAL PROBLEM 


A. CROSS-SECTIONAL HULL 

The mid-section of the imported wireframe hull of the DDG model was utilized 
for this particular study, assuming structural properties remain linear throughout the 
process. Previous work conducted by Newton [Ref. 21] in the area of a two dimensional 
study involved a symmetric cylinder submerged underwater exposed to a plane wave 
incident shock front, utilizing earlier Fortran finite element programs for meshing, 
propagation of underwater shock, and transient response calculations. 

The goal in this particular model is to study responses, specifically peak 
velocities, with and without cavitation effects, of a surface ship model with a draft of 20 
feet. The time frame was limited to the first 30 msecs, in an attempt to eliminate the 
effects of the bubble pulse while comparing plane wave modeling and DAA2 spherical. 

Utilizing the USA code. Fig. 17 illustrates location of cross section in the global 
ship wireframe, 1 foot in depth, and Fig.18 to Fig. 20 illustrate the meshed model. In Fig. 
18, red denotes the structure mesh while green denotes the fluid mesh. Fig. 20 depicts 
the outer DAA, hydrostatic, and structural interface elements mentioned previously. 
Tables 4-6 provide specific model characteristics. 


TOTAL NUMBER OF GRID POINTS 

452 

TOTAL NUMBER OF THIN SHELL ELEMENTS 

230 

THICKNESS OF ELEMENTS 

0.5 in 

TOTAL NUMBER OF WET ELEMENTS 

40 


Table 4. Structure Characteristics 
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MODULUS OF ELASTICITY (ksi) 

29600.0 

DENSITY (Ibm/ft^) 

489.02 

YIELD STRENGTH (ksi) 

51.00 

POISSON’S RATIO(v) 

.3 


Table 5. Element Properties 


NUMBER OF FLUID NODES 3034 

NUMBER OF FLUID VOLUME ELEMENTS 1440 

NUMBER OF DAA FACE ELEMENTS 40 

NUMBER OF HYDROSTATIC PRESSURE FACE ELEMENTS 72 

NUMBER OF PRESCRIBED PRESSURE FACE ELEMENTS 0 
FACE ELEMENTS CONTACTING STRUCTURE 40 

MAXIMUM NUMBER OF FLUID VOLUME NODES 8 

MAXIMUM NUMBER OF FACE ELEMENT NODES 4 


Table 6. Fluid Characteristics 


Utilizing the bulk cavitation program to determine the extent of various bulk 
cavtiation zones of a TNT explosion, the smallest region of bulk cavitation takes place 
directly above the charge at the surface. Therefore, a charge location directly beneath the 
hull at 100 feet standoff (depth = 120 feet) would minimize effects of the bulk cavitation 
zone in the general vicinity of the hull model and provide a means for setting the lower 
edge of the DAA boundary just beyond the cavitating region at a depth of 50 feet, 30 feet 
under the keel. 


B. RESULTS 

Fig. 21 to Fig. 30 record the velocities at the locations in the structure and time 
varying pressures at locations in the fluid volume region for the CFA analysis noted on 
Fig. 19. Fig. 31 to Fig. 33 provide a comparison of keel velocities and accelerations. The 
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pressure profiles, like that seen in the one-dimensional analysis, reflect a respective time 
zero which starts at the point of incident pressure wave arrival at that particular node and 
not the initial time of structure motion. Structure motion is derived from the cumulative 
integration of all the pressure profiles in the fluid region superimposing each pressure 
profile at different spatial locations in time. 

The introduction of a fluid mesh beneath the structure and the turning on CFA 
captured the pressure varying fluid under the hull (Fig. 21 and 24) which significantly 
effected the velocity response of the hull. The cavitation effects on the hull during plane 
wave and DAA2 calculations produced almost identical results in the early time phase of 
hull velocity at the keel (node 42) which was the closest to the charge location. As 
expected in late time response, the plane wave solution begins to slightly vary from the 
DAA2 calculations. However, peak velocity response are the same in early time. 

Pressure profiles at the keel for both CFA models (Fig. 21 and 24), which are 
practically identical, indicate the first signs of cavitation closure pressure occurring at a 
much earlier time compared to the Bleich-Sandler case with an appreciably larger 
pressure peak, 41% of the initial pressure wave. Numerous other cavitation closure 
events occur at the keel with decreasing max pressures. This observed behavior was also 
an event reported in the DDG-53 USS John Paul Jones Live Fire and Testing Shock Post 
report [Ref, 22] which may have contributed significantly to large accelerations at 
particular electrical switchboards following the incident pressure wave. 

Keel velocities in the plane wave (CFA OFF) and DAA2 (CFA OFF) analysis, 
which places the DAA boundary at the ship’s wet hull surface, did not capture the 
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detailed high frequency response of the structure exposed to identical underwater 
explosion conditions (Fig. 27 and 29). Peak keel velocities were identical in magnitude 
and occurrence in time but greatly diverged past reaching the peak velocity of 162 in/sec. 
Pressures along the keel during both CFA OFF calculations were found to include only 
the large initial pressure wave arrival and small variations in the pressure later in time 
along the wet hull, while those in the DAA2 calculation being slightly larger. 

Internal deck velocities, however, revealed velocities greater during the plane 
wave modeling with a surrounding fluid mesh and CFA ON (Fig. 22). The pressure and 
corresponding variations along the wet surface of the hull caused internal velocities to be 
greater than those seen in the DAA2 CFA ON solutions (Fig. 25), primarily due to the 
less accurate calculations late in time for the plane wave equations and higher velocity 
response calculations. Both CFA OFF runs (Fig. 28 and 30) reveal small variations in 
velocity, the greatest occuring at the lowest deck. These small velocites could be 
attributed to the simplified cross sectional hull model which does not include any internal 
longitudinal or transverse bulkheads. If these particular bulkheads were modeled between 
decks, the transit velocity responses of the keel could have been transmitted to the upper 
decks instead of the bending the outer hull only. 

Fig. 31 depicts the overlay of velocities of all 4 models (plane wave with and 
without CFA, DAA2 with and without CFA) at the keel, illustrating the larger velocities 
of the CFA OFF calculations for both plane wave and DAA2 modeling. Plane wave 
modeling of the underwater explosion incident pressure wave without a fluid mesh 
provides for a conservative approach in calculating the peak velocity of the structure at 
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the fluid structure interface. A shift to a higher velocity at an earlier time was more 
pronounced during this particular CFA study in comparison to the Bleich-Sandler case 
with a 22 % increase in peak velocity. Appreciable differences in acceleration were noted 
between CFA ON and CFA OFF calculations (Fig. 32 and 33). CFA OFF results 
reflected initial accelerations 3.5 times larger than CFA ON results but with no following 
variations. 

Peak velocities in each case are provide in Table 7. 


LOCATION 

Node 

Plane wave 
(CFA ON) 

DAA2 
(CFA ON) 

Plane wave 
(CFA OFF) 

DAA2 
(CFA OFF) 

Keel 

42 

124 

124 

162 

162 

Inner Bottom 

3209 

Increasing 

Increasing 

Increasing 

4.3 

Deck 3 

3234 

Increasing 

Increasing 

.094 

.13 

01 Level 

3162 

Increasing 

Increasing 

2.7E-4 

Increasing 


Table 7. Maximum Velocities (in/sec) 


The responses labeled “increasing” in Table 7 reflect increasing values at the end of the 
time period studied. 

Fluid pressures at nodes 3021, 1600,and 2994 during both plane wave and DAA2 
CFA ON calculations (Fig.23 and 26) show similar relative patterns of growing numbers 
and durations of cavitation periods as depth decreases. The magnitude of peak pressmes 
following initial cavitation in the DAA2 results tended to be 15-20% larger than the plane 
wave peak pressures. Node 3021 pressures during both CFA ON calculations at t=.0025 
secs captured the largest peak pressure, followed by Node 2994 at t=.0047 secs. 
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Hence, based on the results, the best conservative approach to modeling a surface 
ship exposed to an underwater explosion would be to apply the plane wave incident 
pressure wave in a fluid mesh that surrounds the entire wet surface of the ship. 
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Figure 18. Structure and Fluid FEM 
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Figure 19. Cross Section Dimensions 
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Figure 20. DAA, Hydrostatic, and Wet Structure Boundaries 
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Figure 21. Velocity and Pressures at the Keel (Node 42)-Plane wave, CFA ON 


52 




















































X 10 


Node 3162 





Figure 22. Internal Deck Velocities-Node 3162(top), Node 3234 (middle), Node 3209 

(bottom).Plane wave, CFA ON 
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Figure 23. Fluid Pressures-Node 3021 (top), Node 1600 (middle), Node 2994 

(bottom). Plane wave, CFA ON 
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Figure 24. Velocity and Pressure at the Keel (Node 42)- DAA2, CFA ON 
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Figure 25. Internal Deck Velocities- Node 3162 (top), Node 3234 (middle), Node 3209 

(bottom). DAA, CFA ON 
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Node 3021 (DAA Boundary) 



Node 1600 (tfid-Fluid) 



Figure 26. Fluid Pressures-Node 3021 (top), Node 1600 (middle), Node 2994 

(bottom), DAA2, CFA ON 
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Figure 27. Velocity and Pressure at the Keel (Node 42)-Plane wave, CFA OFF 
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Figure 28. Internal Velocities-Node 3162 (top), Node 3234 (middle), Node 3209 

(bottom). Plane wave, CFA OFF 
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Figure 29. Velocity and Pressure at the Keel (Node 42)-DAA2, CFA OFF 
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Figure 30. Internal Velocities-Node 3162 (top). Node 3234 (middle), Node 3209 

(bottom). DAA2, CFA OFF 
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Figure 31. Comparison of Keel Velocities. 
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Figure 32. Comparison of Keel Accelerations-Plane wave 
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Figure 33. Comparison of Keel Accelerations- DAA2 
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VI. THREE DIMENSIONAL MODEL 


The ideal three dimensional ship analysis would include: 

1. A detailed ship model which included specific hull and equipment 
characteristics (locations, weights, etc.) 

2. A fluid mesh model of the surrovmding ocean that captures the cavitation region 
in three dimensions, as seen in the one and two dimensional analysis. 

3. A computer hardware that can handle the detailed processing and large solution 
database. 

Due to the complexity of a global ship model and the time integration required, a 
simplified “notional” ship design was utilized which was part of an earlier model of the 
Arleigh Burke Destroyer. Based on the results obtained with the two-dimensional model 
and the requirements for conditions 2 and 3 stated above, a plane wave/CFA OFF 
analysis was utilized for the three-dimensional ship model which sets the DAA boundary 
at the wet surface of the hull. Table. 8 provides general characteristics of the DDG model 
[Ref 23]. 


Length Between Perpendiculars (LBP) 

466.0 ft 

Length Overall (LOA) 

492.1 ft 

Beam, Weather Deck 

66.5 ft 

Draft to Keel 

20.7 ft 

Lightship weight 

6327 LTON 

Full Load displacement 

8000.3 LTON 


Table 8. DDG-51 ASSET Design Summary [Ref. 23] 


Fig. 34 and 35 depict solid modeling of the imported DDG model and specific 
charge locations for the two models created. Fig. 36 illustrates locations of the following 
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major areas of interest were selected for their general purpose of holding vital pieces of 
equipment such as the AEGIS weapon system components and main propulsion units: 
node closest to the charge, bow, stem. Main Engineering Rooms (MER 1 and 2), Combat 
Systems Equipment Rooms (CSER 1, 2, and 3), and the general location of th SPY radar 
room in the superstmcture. 
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Figure 34. Solid Modeling of ASSET DDG-51. 


67 














































Figure 35. Charge Locations Below the Keel (1) and Offset (2). 
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Figure 36. General Space Locations. 
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A. 


1962 FINITE ELEMENT MODEL 


1. Modeling 

This initial model utilized a general plate thickness of 0.4 inches of similar 
material qualities used in the two dimensional model throughout the entire ship. Iimer 
bottom stiffeners and heavy equipment locations were not constructed in the original 
IGES model out of ASSET but were later created in IDEAS for the 5682 element model 
study. Fig. 35 illustrates the 1962 element model which includes 348 wet elements. In 
order to capture peak velocities, damping coefficients were not applied to the model. 

Explosive charge conditions were the same as those utilized for the two 
dimensional analysis: 


CHARGE TYPE 

TNT 

WEIGHT 

100 LBS 

LOCATION (relative to ship) 

MIDSHIPS, CENTERLINE, 100 ft beneath the keel 

DEPTH 

120 feet 


Table 9. Charge characteristics 


2. Results 

The whipping motion of this particular model exhibited the expected basic motion 
of a simple ffee-ffee beam excited by a forcing function in the middle of the model. The 
highest velocities and accelerations occurred at the node (Node 42) closest to the charge 
(Fig. 37 and 38) due to an unsmeared and unstiffened keel and outer hull. The 
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unsmeared and unstiffened keel also experienced high frequency oscillations in 
comparison to other areas of the model. The bow and stem experienced the second 
largest maximum vertical velocities and accelerations (Fig. 39 to 42) which occurred 
approximately 45 msec and 32 msec, respectively, after the peak keel velocity. 

Peak velocities at the bow,stem, and closest point to the charge at the keel indicate 
the greatest magnitude in velocities to be in the vertical direction, followed by fore and 
aft motion, and the athwartships direction to be the lowest. However, internal deck 
locations exhibited the greatest peak velocities in the fore and aft direction, followed by 
athwartship, and vertical velocities to be the lowest because of their general locations 
close to the incident shock wave arriving midships in the model. Both MER spaces (Fig. 
43 to 46) experienced the largest internal accelerations, specifically in the fore and aft 
direction followed by CSER 3 (Appendix C) with almost equal velocites and 
accelerations in the fore and aft and athwartship directions. MER 2 experienced the 
largest fore and aft motions in comparison to the other internal spaces studied, largely due 
to its close location to the initial shock wave. Amongst the combat systems spaces, the 
SPY radar room (Appendix C) had the highest fore and velocity (9 in/sec) but with 
similar accelerations. The asymmetry of the aft super structure arrangement places a 
majority of the stmcture to port. This asymmetric condition in the DDG-51 model 
contributed to the athwartship velocities and accelerations in all areas of the ship. 

Modeling the ship globally without a fluid volume mesh and exposing it to a 
plane wave vmderwater explosion captured the high frequency vibration or “ringing” in 


71 


the hull which was not observed in the plane wave/CFA OFF two dimension^ cross- 
sectional hull model. Table 10 summarizes maximum velocities and accelerations: 


Location 

Node 

Fore and Aft 

Athwartships 

Vertical 



V(in/sec) 

A(G) 

V(in/sec) 

A(G) 

V(in/sec) 

A(G) 

Keel 

42 

10.2 

47.1 

3.3 

20.8 

176 

2210 

Bow 

842 

13.1 

9.2 

9.1 

. 4.3 

27.2 

10.1 

Stem 

2157 

13.4 

10.1 

3.8 

1.33 

20.1 

12.2 

MERl 

1933 

7.2 

12.1 

2.6 

2.2 

.7 

.058 

MER2 

1895 

12.3 

9.92 

4.6 

3.8 

4.2 

.13 

CSERl 

1649 

5.5 

3.6 

2.5 

1.3 

.17 

.019 

CSER2 

1748 

6.5 

5.1 

1.33 

1.1 

.16 

.025 

CSER3 

1990 

6.1 

7.1 

6.8 

5.03 

.135 

.08 

SPY Radar 
Room 

2534 

9 

4.8 

2.6 

1.8 

.055 

.0088 


Table 10. Maximum Absolute Velocity and Acceleration (1962 Elements) 
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Figure 37. Keel Velocity (Node 43) 
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Figure 38. Keel Acceleration (Node 43) 


74 























































Velocity (in/sec) Velocity (in/sec) Velocity (in/sec) 


Fore and Aft Direction 





Time(sec) 


Figure 39. Bow Velocities (Node 842) 
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Figure 40. Bow Accelerations (Node 842) 
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Figure 41. Stern Velocities (Node 2157) 
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Figure 42. Stern Accelerations (Node 2157) 
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Figure 43. MER 1 Velocities (Node 1933) 
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Figure 44. MER 1 Accelerations (Node 1933) 
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Figure 45. MER 2 Velocities (Node 1895) 
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Figure 46. MER 2 Accelerations (Node 1895) 
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B. 5682 FINITE ELEMENT MODEL 


1. Modeling 

In order to provide a detailed mesh for the wet surface of the hull and improve 
equipment locations, a finer mesh was created of the same model which included inner 
bottom stiffeners and 955 wet elements, Fig. 47. Utilizing the available ASSET data on 
main engineering equipment [Ref 23], early Flight II/IIA study drawings [Ref 24], and 
DDG combat system equipment data [Ref 25], large groups of heavy equipment were 
lump massed in their respective areas. 

The detailed requirement for longitudinal and transverse stif&iers was simplified 
by the smearing technique utilized in ASSET (APPENDIX D). Although previous work 
by Cunningham [Ref. 26] revealed the limitations in local and global smearing, the 
ASSET technique sizes stiffeners to provide a T-section with reasonable proportions that 
meet strength and geometric constraints. The stiffener cross-sectional area is smeared into 
the plate cross sectional area, creating a smeared plate thick than the standard plate by a 
factor of (1 + smear ratio(S)). Based on the data for the notional DDG design which 
included intitial plate thickness and stiffener smearing ratios, the following plate 
thicknesss were used in this particular model: 
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Structure Location 

Smeared Plate Thickness (in) 

BOTTOM SHELL 

.57 

INNER BOTTOM 

.42 

DECKl 

.51 

DECK 2 

.54 

DECK 3 

.49 

DECK 4 

.48 

MAIN DECK 

.55 

TRANSVERSE BULKHEADS 

.33 

OUTER HULL 

.47 

SUPERSTRUCTURE 

.33 


Table 11. Smeared plate thicknesses (inches) 


In addition to Appendix E which denotes general weights and locations of 
combat system equipment that were lump massed, the following specific engineering 
room equipment was used with respect to their locations derived from [Ref. 24]. 


GAS TURBINES (4) 

84.7 

GAS GENERATORS (3) 

94.2 

REDUCTION GEARS (2) 

109.5 


Table 12. Engineering Equipment Total Weights (LTON)[Ref. 23] 


The net weight of the entire model was 3598 Eton, 43% lighter than the original 
“light weight” DDG ASSET model. 

Charge location was placed in an offset position utilizing only 120 lbs TNT at a 
standoff distance of 107 feet and a depth of 50 feet . 
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Figure 47. DDG 5682 Finite Element Model 
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2 . 


Results 


Due to the increased complexity of the meshed model, post processing of this 
particular model, became quite difficult during the MSC/NASTRAN restart which 
prevented further post-processing into IDEAS. Fig. 48 to 51 depict the greatest velocities 
and accelerations occuring midships of this model, specifically at the point along the hull 
closest to the charge (Node 393) and at the node along the keel with the same axial 
distance from the bow (Node 397). These two locations coincide with those same points 
on the hull noted previously in Fig. 35 with charge 1 for node 397 and charge 2 for node 
393. 

Bow and stem reactions (Fig. 52 to 55) were noted higher than those seen in the 
1962 element model due to the majority of concentrated masses being located in central 
locations of the ship (i.e. aft of the forward VLS and forward of CSER3) as well as the 
spaces of interest (Fig. 36) and the lack of masses in the general area of the bow and 
stem. Peak bow and stem vertical reactions occurred 40 msec and 30 msec, respectively, 
later relative to midships reactions (Node 393 and 397) with high frequency oscillations. 
MER 1 (Fig. 56 and 57) exhibited lower velocities and accelerations relative to MER2 
(Fig. 58 and 59) which can be contributed by the asymmetry between the two 
engineering rooms. The gas turbine modules and reduction gear masses were located on 
the port side of the engine room of MER 2 which was the same side as the charge 
location and, as seen in the 1962 element model, MER 2’s location was closer to the 
incident shock wave arrival point along the hull. Bow,stem, and both MER spaces tended 
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to have reactions in the following order of highest to lowest magnitude: vertical, fore/aft, 
and athwartships. Combat system equipment rooms (CSER 1,2 and 3) and the SPY radar 
room (Appendix F) showed very little motion in the fore/aft and athwartships direction 
with practically zero vertical motion. 

For this particular tactical scenario and model, the only shipboard internal space 
which experienced the highest acceleration was MER 2 with large magnitudes of 
oscillations in the fore and aft direction as well as athwartship (Fig. 58 and 59). All other 
internal spaces encountered low accelerations which did not seriously introduce severe 
shock loadings upon their respective equipment. The effect of smearing and lump 
massing groups of large pieces of internal equipment created dramatically different 
results compared to the 1962 unsmeared, unlumped element model. In addition to these 
structural differences, the low charge weight and offset distance of the underwater 
explosion did not produce appreciable effects on internal ship behavior. Table 12 
provides a summary of the maximum velocities and acclerations in each area: 
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Location 

Node 

Fore and Aft 

Athwartships 

Vertical 



V(in/s) 

A(G) 

V(in/s) 

A(G) 

V(in/s) 

A(G) 

Standoff 

Point 

393 

17.1 

64.1 

53.2 

172 

58.9 

210 

Keel 

397 

10 

52 

30.2 

no 

28.8 

106 

Bow 

1554 

13.2 

42.3 

12.1 

44 

21.2 

58 

Stem 

2281 

13.8 

23.1 

10.8 

30.1 

12.2 

35.6 

MER 1 

4448 

.33 

.18 

2.1 

.67 

2.23 

1.08 

MER2 

4543 

8.8 

42.4 

9.2 

44.5 

39.3 

133 

CSERl 

3567 

.013 

1.9E-3 

.012 

1.3E-3 

6.4E-4 

5.1E-5 

CSER2 

3669 

.19 

.023 

.78 

.07 

7.7E-4 

1.48E-4 

CSER3 

3541 

.15 

.19 

1.4 

.23 

7E-3 

1.5E-3 

SPY Radar 
Room 

6955 

.59 

.071 

.15 

.019 

3.2E-4 

4.3E-3 


Table 13. Maximum A 


3solute Velocities and Accelerations (5698 element model) 
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Figure 48. Hull Standoff Point Velocities (Node 393) 
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Figure 49. Hull Standoff Point Accelerations (Node 393) 
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Figure 50. Keel Velocities (Node 397) 
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Figure 51. Keel Accelerations (Node 397) 
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Figure 52. Bow Velocities (Node 1554) 
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Figure 53. Bow Accelerations (Node 1554) 
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Figure 54. Stern Velocities (Node 2281) 
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Figure 55. Stern Accelerations (Node 2281) 
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Figure 56. MER 1 Velocities (Node 4448) 
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Figure 57. MER 1 Accelerations (Node 4448) 
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Figure 58. MER 2 Velocities (Node 4543) 
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Figure 59. MER 2 Accelerations (Node 4543) 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The introduction of a fluid volume model with a structure and utilizing the USA 
/CFA interface provides a means to capture the effects of cavitation upon a surface vessel 
in the early time frames of an underwater explosion. The effects of cavitation closure as 
seen in the two dimensional model can be quite significant and inflict another impulsive 
force upon the ship as reported in the DDG 53 shock trial reports [Ref. 22] in addition to 
the incident wave, bubble pulsations, bottom reflection, or bottom refraction. This 
increases the complicated task of predicting surface ship behavior when exposed to an 
underwater explosion. CFA OFF modeling may provide a conservative approach in 
predicting ship response to an imderwater explosion in early time. Finite element 
modeling, time integration calculations, and tactical placement of an explosive device 
must be carefully done to accurately model surface ship behavior. Plane wave modeling 
with a surface ship model and a fluid mesh may provide more accurate overall ship and 
internal equipment responses. 

With improving systems and capacity, the testing of new ship models with 
specific interior designs from the ship designer with the underwater shock analyst 
utilizing the USA/CFA code can prove to be a formidable method in predicting forcing 
fimctions on internal equipment of a surface vessel. In view of the changes of the DDG 
Flight IIA designed compared to the present DDG Flight II, this particular method may 
provide some cost benefits in future ship shock trials. 
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Recommend continuing studies in the following areas: 

1. Improve on the two and three dimensional model with internal masses or fluids. 

2. Obtain actual DDG-51 Flight IIA DDG designs from NAVSEA which can be 
imported into an FEM program such as IDEAS or MSC/PATRAN which can be 
tested with the USA code for both linear and non-linear responses. 

3. Introduce other surface ship model types such as the slice hull. 

4. Model the underwater tactical scenario with significant shallow water bottoms 
in order to study “littoral” responses of a surface ship close to the shore. 
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APPENDIX A. BULK CAVITATION PROGRAM 


The following program calculates the bulk cavitation zone for various weights and 
depths of TNT. The program can be easily modified for other explosive material. 
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%The following program will compute the bulk cavitation zone of an underwater 
%explosion of TNT explosive with depth and weight of charge as inputs. 

% 

%TNT characteristics 
A1=1.18;A2=-.185;K1=22505;K2=.058; 


%D=input('Enter depth of your charge ') 

%W=input('Enter the weight of TNT explosive ') 
Co=4967; 

Pa=14.77 

x=(0:5:l200] ' ; 
n=length(x); 
yu=diag(zeros (n)) ; 
yl=diag(zeros(n)); 
for j=l:n 

yd) =0; 

for i=l;500 

rl = sqrt ( ( (D-y ( i) ) ^^2) ^ (x { j ) ^^2 ) ) ; 
r2=sqrt( ( (D+y(i))''2) + (x(j)''2)); 
theta=K2*W^(l/3) M (W'' (1 /3 ) /r 1) ^A2 ) ; 
dr= (r2-rl)/(Co*theta/1000) ; 

Pmax=Kl* ( {W^(l/3) /rD'^Al) ; 

Pl=Pmax*exp(-dr); 

Pl = Pinax*exp (-dr) ; 
P2=KlM(W''(l/3)/r2)''Al) ; 

Ph=64*y(i)/144; 


%speed of sound in water (ft/sec) 
%atmospheric pressure 


%real charge distance 
%image charge distance 
%decay constant (msec) 

%decay rate 
%Incident shock wave 

%Inc shock wave @ raref arriv 
%Inc shock wave eraref arr 
%rarefrac at arrival time 
%hydrostatic pressure 


F = Pl-^Pa+Ph-P2 ; 

if (F<=0) 

break 

else 

y ( i + 1)=y(i)+.1; 

end 

end 

yu (j)=y(i); 
y{l)=0; 

for s=l:500 

rl = sqrt ( { (D-y (s) )''2) + (x(j)'‘2) 
r2=sqrt { ( (D-^y (s) )"2) + (x(j)''2)) ; 
theta=K2*W'' (1/3) * ( ((1 / 3 ) / r 1) A2 ) ; 
Pmax=Kl* ((W''(l/3)/rl)''Al); 

Pl=Pmax*exp(-dr); 

P2=K1*((W^(1/3)/r2)^Al); 

Ph=64*y(s)/144; 

dr=(r2-rl)/(Co*theta/1000); 

gl=(D^y(s))/r2; 

g2=r2-(2*D*gl); 

g3={g2/rl)*{((A2*r2)/rl)-A2-1); 
g4=Pl/(Co*theta/1000); 
g5= {A1*P1) / (rl''2) ; 
g6=Pl+Pa+Ph; 

G=- (g4*(l+g3))-(g5*g2) + {64*gl/144) + 
if {G>=0) 
break 
else 

y (s + l) =y (s ) +•. 1 ; 

end 

end 

y1(j)=y (s) ; 

if (yl{j)<=yu(j)) 

z = j; 

break 

else 

2 = n ; 

end 


%F{x,y) equation for upper boundary 


%real charge distance 
%image charge distance 
%decay constant (msec) 

% Incident shock wave 
%Inc shock wave @ raref arriv 
%rarefrac at arrival time 
%hydrostatic pressure 
%decay rate 


(Al*g6/r2); %Combined G{x,y) lowcav bound 


end 

yu=('l)*yu(l:z); 
yl=(-l) ‘yl(l:z) ; 
xf=x{1;z); 
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APPENDIX B. CROSS-SECTIONAL HULL INPUTS 


The following are sample inputs for the two-dimensional model only. 


assign outpuc^^' UTi ' , unir:=li, aeiece i? 
assign output2='IDEASl.op2",unit=12,unknown 
ID USA,CROSS HULL3 
SOL 109 $ MSC 

DIAG 8,13 $ 

TIME 99 
$ 

$ START SDRC'S I-DBAS OPTION 

compile phaseO souin=mscsou nolist noref $ 

$ add CALL to get output2 geometry datablocks early in solution, 
alter 297 $ before call seprep2 
type parm,,i, y,ounit2 = 12 $ 
if (post < 0) then $ 

OUNITl=OUNIT2 

CALL OUT2GEOM CSTM, GPL, GPDT, EPT, MPT, GEOM2 , GEOM3 , GEOM4 , 
GE0M1,DIT, DYNAMICS, EQEXIN, BGPDT, CASECC// 
OUNITl/0 $ 

endif $ 

$ END SDRC'S I-DBAS OPTION 
$ 

compile subdmap=selg,souin=mscsou,noref,nolist $ 

ALTER 14 $ AFTER SSGl 

type parm,,i,n,storit = 0 

call dbstore pjx,,,,//seid/O/'DBALLVs,storit $ 

OUTPUT2 BGPDTs,EQEXINs,,,///ll//25000 $ 

$ 

compile subdmap=sekr,souin=mscsou,noref,nolist $ 

ALTER 49 $ AFTER gpsp 

type parm,,i,n,foundit = 0 $ 

call dbfetch /pjx,,,,/seid/O/O/O/s,foundit $ 

OUTPUT2 uset,pjx,cstms,,///ll//25000 $ 

$ 

compile subdmap=DTRANRS,souin=mscsou,noref,nolist $ 
alter 1 
type db ECTS 

ALTER 6 $ before TRDl 

OUTPUT2 ECTs,BDD,KDD,MDD,///11//25000 $ 

$add EOF marker 
OUTPUT2 //-9/12/ $ 

$ 

EXIT $ 

CEND 

TITLE=CROSS HULL 

LOADSET=10 

TSTEP=25 

SPC = 1 

DlSP(plot)=ALL 

velocity(plot)=all 

stress(plot)=all 


BEGIN BULK 

PARAM AUTOSPC YES 

PARAM POST -2 


$ TSTEP,SID,N1,DT1,N01 
TSTEP,25,40,0.025,1 
$ 

$ LSEQ,SID,DAREA,LID 
LSEO,10,10,1 

GRID 1 0 -336-0 

GRID 2 0 -336.0 


240.0 12.0 0 

240.0 0.0 0 
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CROSS SECTIONAL HULL3B PROBLEM 

1 , 0 , 0 , 1 , 1 , 0,0 


USAF 

GAL 

F F F F 
F F F 


AUGMAT INPUT FOR CROSS HULL3 PROBLEM 
GAL HULL3.FLU HULL3.GEO HULL3.PRE 


$ STRNAM FLUNAM GEONAM PRENAM 


F F F F 

$ 

FRWTGE 

FRWTST 

FRWTFL 

LUMPFM 

F F F T 

$ 

FLUSKY 

DAAFRM 

SYMCON 

DOFTAB 

T T F F 

$ 

PRTGMT 

PRTTRN 

PRTSTF 

PRTAUG 

F F F F 

$ 

MODTRN 

STRLCL 

INTWAT 

CFADYN 

NASTRAN-MSC 

$ 

MAINKY 




0 

$ 

NTYPDA 




452 1356 6 3 

$ 

NSTR 

NSFR 

NFRE 

NFTR 

1 

$ 

NSETLC 




0 1 80 1 

$ 

NDICOS 

JSTART 

JSTOP 

JINC 


FLUID MASS RUN FOR CROSS SECTION HULL3B 
HULL3.FLU HULL3.GEO GAL HULL3.DAA 
T 


FLUID.VOL FLUID.VOL 
NASTRAN-MSC 
0 452 0 80 
0 0 40 
0 3034 

9.3455-5 57120. 

10 

240 . 0. 1. 0. 

14.7 386.4 
0 
0 
0 


$ FLUNAM GEONAM GRDNAM DAANAM 
$ PRTGMT PRTTRN PRTAMF CALCAM 
$ EIGMAF TWODIM HAFMOD QUAMOD 
NASTAM STOMAS STOINV 
FRWTGE FRWTGR FRESUR 
STOGMT ROTGEO ROTQUA 
STRMAS SPHERE ROTSYM 
CAVFLU FRWTFV INTCAV 
MASREF 
FACNAM 


PCHCDS 
FRWTFL 
RENUMB 
PRTCOE 
OCTMOD 
BOTREF 
CAVNAM 
MAINKY 
$ NSTRC NSTRF NGEN 
$ NBRA NCYL NCAV 

$ NFVNC NFVNF 

$ RHO CEE 

$ NVEC 

$ DEPTH CXFS CYFS 
$ PATM GRAVAC 

$ NUMCON 
$ NSRADI 
$ NSORDR 


NGENF 


CZFS 
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USA-NASTRAN CFA TIMINT RUN FOR CROSS HULLS PROBLEM 


HULL3.PRE HULL3 

. POS 

$ PRENAM POSNAM 


HULL3.RES 


$ RESNAM WRTNAM 


F F F 

$ 

REFSEC 

FLUMEM 

STRVEL 


1 0 

$ 

NT I NT 

NCHGAL 



0. 3.E-5 

$ 

STRTIM 

DELTIM 



3.E-2 0. 0. 

$ 

FINTM 

ADAMP 

BDAMP 


T F F F 

$ 

EXPWAV 

SPLINE 

VARLIN 

PACKET 

F T F F 

$ 

HYPERB 

EXPLOS 

DOUBDC 

VELINP 

F F F 

$ 

BUBPUL 

SHKBUB 

BOTREF 


1 

$ 

NCHARG 




0. 14.7 

$ 

HYDPRE 

PATM 



0. -1.44E3 0. 

$ 

XC 

YC 

ZC 


0. -351.8 0. 


$ SX 

SY 

sz 


201 

$' 

JPHIST 




1. 0. 

$ 

PNORM 

DETIM 



3.E-5 

$ DTHIST 




2 

$ 

CHGTYP 




100. 100. 120. 


$ WEIGHT SLANT CHGDEP 

1 150 0 

$ 

NSAVER 

NRESET 

NSODFL 


0 0 0 0 

$ 

LOCBEG 

LOCRES 

LOCWRT 

NSTART 

T F F 

$ 

FORWRT 

STBDA2 

ASCWRT 


T F F 

$ 

NOAMBI 

PRTVOL 

PRTINI 


o 

in 

$ 

FVBETO 

FVBETl 



111 

$ 

ICAVSW 

lORDF 

lORDV 


124 1440 3034 

$ 

NFVWAV 

NFVELM 

NFVNOD 


T 

$ 

DISPLA 




0 0 

$ 

NPREVT 

NPREVF 



T T 

$ 

LISTRE 

PRTPLT 



2 0 0 

$ NWETHS NDRYHS NUMSET 


42 2 

$ 

NODOUT 

NFROUT 



3209 2 

$ 

NODOUT 

NFROUT 



2 0 0 

$ 

NWETHS 

NDRYHS 

NUMSET 


42 2 

$ 

NODOUT 

NFROUT 



3209 2 


$ NODOUT NFROUT 


4 0 

$ 

NPREHS 

NUMSET 



42 0 

$ 

NEQHPR 

IPRS 



1600 1 

$ 

NEQHPR 

IPRS 



2994 0 

$ 

NEQHPR 

IPRS 



3021 0 

$ 

NEQHPR 

IPRS 



F 

$ 

SCALEF 





hullSB MODEL 
HULLS.POS 
1001,40,0,0,1 
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assign inputt2='UTl',unit=ll $ 
assign inaster=' hullSB .MASTER' $ 
assign output2='IDEAS2.op2',unit=12,unknown 
RESTART version=l/ keep $ 

ID USA, HULL3B 
SOL 109 $ MSC 

DIAG 8,13 $ 

TIME 120 

$ 

compile subdmap=DTRANRS,souin=mscsou,noref $ 
ALTER 7,7 $ eliminate TRDl 
delete /PNLD,,,,/ $ 

INPUTT2 /UHV,,,,//ll $ 

EQUIVX UHV/UDVT/-1 $ 

CEND 

TITLE= THESIS HULL3B MODEL 

LOADSET=10 
TSTEP=46 
DISP(plot)=ALL 
STRESS(plot)=ALL 
VELOCITY(plot)=ALL 
ACCELERATION(plot)=A11 
BEGIN BULK 

TSTEP,46,1001,3.Oe-5,1 
ENDDATA 
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APPENDIX C. COMBAT SYSTEMS AREA RESPONSES-1962 ELEMENT 

MODEL 


The following plots represent the responses for CSER 1,2,3 and the SPY radar 
room for the 1962 element model. 


/■ 


109 


Velocity (in/sec) Velocity (in/sec) 


Fore and Aft Direction 




Figure 60. CSER 1 Velocities (Node 1649) 
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Figure 61. CSER 1 Accelerations (Node 1649) 
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Figure 62. CSER 2 Velocities (Node 1748) 
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Figure 63. CSER 2 Accelerations (Node 1748) 
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Figure 64. CSER 3 Velocities (Node 1990) 
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Figure 65. CSER 3 Accelerations (Node 1990) 
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Figure 66. SPY Radar Room Velocities (Node 2534) 
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Figure 67. SPY Radar Room Accelerations (Node 2534) 
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APPENDIX D. SMEAR PLATE DATA 
The following illustrates the smeared plate area method and ratios that were 


extracted from [Ref. 23], 


SMEARED PLATE 


STIFFENED PLATE 


Aj (Area of Tee) 



f 


Figure 68. Smeared Plate Area Method [Ref. 14] 


SEGMENT PROPERTIES 


■AREA 



TOTAL 

SHEAR 

SEG 

IN2 

IN2 

1 

11.00 

.66 

2 

8.13 

. 72 

3 

8.01 

. 71 

4 

10.01 

1.06 


PROPERTIES OF 
N.A. TO 
PLATE 
IN 

•51 28.56 

•69 19.71 

•68 19.78 

•95 36.86 


SMEAR 

FLANGE WT/FT RATIO 

IN3 LBF/FT 

2.96 37.37 .10 

3.86 27.61 .22 

3.84 27.22 .22 

6.64 33.99 .22 


STIFFENED PLATES 

-SEC MOD- 

PLATE 
IN3 
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APPENDIX E. COMBAT SYSTEM EQUIPMENT WEIGHTS 


The following table comprises the general combat system equipment utilized in 
the distribution of lump massing in the 5698 element model whose weights derived from 
[Ref. 25] and placed in locations utilizing [Ref 24]. 


Space 

Location 

Weight (lb) 

Communication Center 

2-126-1-C 

1200 

Radio Trans. Room 

2-158-1-C 

5884 

Array Room No. 1 

03-128-2-Q 

14692 

Array Room No. 2 

03-128-2-Q 

14692 

Array Room No. 3 

03-155-1-C 

14692 

Array Room No. 4 

03-155-2-C 

14692 

Radar Room No. 1 

03-128-0-C 

27720 

Radar Room No. 2 

03-142-0-C 

27720 

Radar Room No. 3 

01-274-1-C 

12035 

Power Supp/Conv Room 

3-126-0-Q 

11850 

Power Conversion Room 

3-319-0-Q 

11850 

Sonar Equip. Room No. 1 

1-18-0-Q 

11545 

Sonar Equip. Room No. 2 

2-18-0-Q 

11545 

Sonar Equip. Room No. 3 

3-18-0-Q 

11545 

Sonar Dome 

External to bow 

117246 

Combat Sys. Equp. Rm 1 

2-53-1-C 

4000 

Combat Sys. Equp. Rm 2 

2-126-2-C 

23285 

Combat Sys. Equp. Rm 3 

1-300-0-C 

4000 

Combat Inform. Center 

1-126-0-C 

12120 

Tomahawk Equp. Room 

2-153-2-C 

1421 

5”54 Loader Drum Room 

1-46-0-C 

41658 

VLS No. 1 

01-78-0-M 

138453 

VLS No. 2 

01-338-0-M 

258853 

Starboard Torpedo Tube 

01-342-1-M 

2430 

Port Torpedo Tube 

01-342-2-M 

2430 

Gun Mount 21 

03-115-0-M 

13600 

Gun Mount 22 

02-310-0-M 

13600 
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APPENDIX F. COMBAT SYSTEMS AREA RESPONSES-5698 ELEMENT 

MODEL 


The following plots represent the responses for CSER 1,2, 3 and the SPY radar 
room for the 5698 element model. 
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Figure 69. CSER 1 Velocities (Node 3567) 
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Figure 70. CSER 1 Accelerations (Node 3567) 
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Figure 71. CSER 2 Velocities (Node 3669) 
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Figure 72. CSER 2 Accelerations (Node 3669) 
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Figure 76. SPY Radar Room Accelerations (Node 6955) 
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